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LIST OF SYMBOLS 


Capacitors for the storage of the desired heat flow signal 
Capacitor for the storage of the radiant heat flow signal 


Capacitor for the. storage of the adiabatic wall temperature 
signal 


Control voltage analogous to (qd, - a) 
Coefficient of heat transfer (BTU/nr /ft” /°F) 
Stefan Constant (pru/nr /et* /°F abs) 


Thermal conductivity (Bru/ft"/°F) 


Amplifier gain 
Convective heat flow (BTU/hr /ft") 
Radiant heat flow (Bru/nr /ft") 


q + qd. = desired heat flow (pru/nr /£t”) 


Resistor (ohms) 


Feedback Resistor 
Adiabatic wall temperature (°F abs) 
Surface temperature (°F abs) 


Time 


High gain operational amplifier 
Operational amplifier used as a peak follower 
Operational amplifier with feedback 


Output voltage 


Input voltage 


- Y= 





I. SUMMARY AND INTRODUCTION 


During the past decade aircraft speeds have extended further and 
further into the supersonic regime’. Today, most aircraft design 
companies project civil transports to fly between Mach 2 and Mach 3. 
Such machines will experience structural problems due to the aerodynamic 
heating associated with these high speeds. Therefore there is a need 
for structural research into the effects of this type of environment. 
Early studies on radiant heat techniques for thermal simulation were 
made by Abramson in the United States and Horton! 2222524351, ,, Great 
Britain concurrently. Their works proposed computer controlled systems 
using quartz tubular heaters as the basic heat producing elements. The 
concept of a time-shared computer for the calculation of desired radiant 


[5], 


heat flow was described by Horton The idea has been amplified and 


enlarged in this work. 
A practical aerodynamic heating computer-controller must be: 


(1) Capable of producing the correct temperatures over a given 


range of Mach numbers. 
(2) Relatively simple to operate. 
(3) Dependable. 
(4) Safe to operate. 
_(5) As inexpensive as possible. 
To create the thermal environment corresponding to a given Mach 


number, the simulator must produce the appropriate heat flow. This has 


three components. 
(1) The convective heat transfer to the body. 
(2) The radiant heat transfer to or from the body. 


(3) The internal heat generated within the body. 


é 





Xe . 
Superscript numerals in brackets refer to the References at the 
end of the report. 





@ @ ’ 


It is normally accepted that the Newtonian law can be used to compute 
the convective heat Flow. 19! Thus, the convective heat input condition 


is defined by 


q=h(t. - t,) 
where :' 
: 2 
q = heat flow per unit area (BTU /ft ) 
h = heat transfer coefficient (ptu/nr/et* /2F) 
ay = adiabatic wall temperature (°F abs) 


te = surface temperature (ant abs ) 


The purpose of this work was to design and build a control system 
to meet the above conditions. Even though the equipment is the result 
of a university research project, it is felt that many of the ideas 


Bacveloped are directly applicable to industrial plants. 


The results of the investigation show that the use of a time-shared, 
digital-analog computer control is feasible and most desirable from a 
standpoint of cost, accuracy, simplicity and dependability. Although 
the control was limited to a four channel heating system, the results 
show that it is possible to expand the system to more channels, and to 


include loading should this be desired. 





II. BASIC LOGIC AND BLOCK DIAGRAM 


As stated previously the accepted Newtonian equation: 


qa = h(t. : t.) 


was used to determine the desired convective heat flow, q, into the 
body. External radiant heat input, 4q_ (such as from a bomb flash or 
solar radiation) was considered to be a quantity which could be positively 
defined and supplied as an analog from a function generator. (In fact, 
the basic digital-to analog converter was adapted for this: purpose. ) 
Radiant heat loss in flight was considered as a correction term to the 


q supplied, since the computation K t, and the computation 


i 


Ot . 
Ke ot 


are closely allied. 


The quantities ba h , and q. Can, for*any specifics lighy 
program, be defined explicitly as time dependent variables. Therefore 
analogous time dependent voltages can be produced from appropriate func- - 
tion generators. However, it is our premise that if the operations 
involving Cae » hy, and q. can be sequentially programmed, then a 
Single digital to analog converter can be adapted to generate all. 

(See Fig. 1) In the system developed, a tape reader ‘is used to set up 

a voltage analogous to oe at time t in the program. This voltage is 
stored in a simple capacitor store using a peak follower circuit as 
described later. Next, at the time t + At the adiabatic wall temper- 
ature, eee » is generated. The appropriate surface temperature, ee 
is measured by a Piernocouciend and is combined with Cae in a summa- 
tion circuit described in Section III. The output from the summer is 
stored and when the tape advances to read the heat transfer coefficient, 
the stored temperature potential is fed into the converter as a reference 
voltage. The converter then multiplies the tempefature potential by the 
heat transfer coefficient, thus computing aq. At the same time ae 
previously stored, is added to q giving a voltage analogous to the 


desired heat flow, » uncorrected for radiant loss. The signal is 


qa 





fed into a comparator where it is compared with the measured heat flow 
corrected for radiant loss. The difference, C , between the computed 
and measured values is used to adjust the output of the saturable core 
reactors so that C converges to the control point for correct heat 
flow. Details of the comparison system, the saturable core reactor 
control, and the heat input determination are given in Ref. [8] and [9] 


respectively. 


The Nacte logic shown in the block diagram (Fig- 1) was used in the 
design of the control system discussed. In order to satisfy the design 
requirements as previously stated, the system incorporates time sharing 
and digital to analog conversion techniques. Figure 2 is a diagram of 
the time-shared computer control system. In order to familiarize the 
reader with its operation, the time-wise sequence of events required to 
compute qq? will be followed. 


(The system components are described in Section III, and Appendix 


A contains basic analog computer Pecrmicnees Lo 


First, assume the wipers are on pole one of each of the nine banks 
and the tape input command is a The reader generates a voltage 
proportional to ca from the reference voltage and this is fed through 
pole one of bank one to the storage capacitor si - The switch is now 
stepped to position two and the tape command is simultaneously changed 
vie) ces - The reader generates a voltage analogous to coe from the 
original meLerence: VOlLAase-™ L » aS a voltage, passes through pole 


aw 


two of bank one to the storage capacitor C The output of peak 


follower two is analogous to eee and, Es four, has we sub- 
tracted. This quantity, eae ~ t)s is now an analog signal appearing 
at pole three of bank one. Once again the switch advances and the tape 
command changes to h. The analogous voltage to tee - t.) appears 
as a reference voltage and is multiplied by the command signal in the 
reader head giving q. At the same time, q. is added to q in ampli- 
fier one giving the analogous voltage of qa.” This voltage passes 
through pole three of bank six, seven, and eights and is stored in Coa 
and C, - The output of peak follower six is a continuous voltage 
analogous to dai ° As the steps continue, channels two, three, and 


four are filled in order. 





Banks four and five are used to empty the storage capacitors in 


order to allow the storage of decreasing voltages. 


Banks seven and eight isolate Cy - this is necessary in order 


‘to prevent discontinuous signals of ds and still allow for decreasing 


voltage levels analogous to desired heat flow. 





III. COMPONENT DESIGN EVALUATION 


An eight pin Soroban tape reader is utilized in the function gen- 
erator section of the control system. The reader is used to connect 
various combinations of resistors to amplifier one. As shown in 


Appendix A, the output voltage of a summer amplifier is: 


Hale: R R 
or (at te += + a] Ve 
ioe Ss 


Therefore, by proper choice of resistor values, and by correct choice 

of hole combinations in the tape, the reader head in conjunction with 
amplifier one, will produce any multiple of a given voltage. The prin- 
ciple of this operation is illustrated in Fig. 3, a simplified diagram 
of the reader head amplifier combination. The choice of the resistors 
was made to correspond with the values of numbers on a Soroban tape 
puncher used in programming. Since the tape puncher keyboard is limited 


to numbers from 1 to 159, the maximum inaccuracy of generation is 


Xx 


~ B(159) 


. 


where X is the operating reference voltage range. If full utilization 
of the eight holes had been made in the tape puncher, 255 combinations 


would result, with a maximum inaccuracy of 


va 


Xx 


* B(255) 


Therefore, unless conformity to available equipment is a criteria, any 


accuracy desired may be obtained, since the inaccuracy is equal to 


+ x 
Pe en 
where n is the number of holes. 


Amplifiers used in the digital-analog computer control are Model 
K2-XA operational amplifiers made by Philbrick Researches Incorporated. 
The K2-XA has a small Signal rise time of 1 Usec , an input impedance 


26 





of over 100 Meg 2, output range of #100 v , and a drift rate referred 
to the input of approximately 8 mv per day. 


Amplifier 3, Fig. 2, is a chopper stabilized amplifier with a DC 
gain of 1000. This amplifier uses a Philbrick K2-P stabilizing amplifier 
in conjunction with a K2-XA in order to increase stability, giving it a 


long term drift rate of Jess than 0.1 mv /day referred to the input. 


The peak follower circuit mentioned previously is simply a K2-xXA 
amplifier with a feedback loop of zero resistance. <A capacitor, isolated 
from the rest of the circuit by a diode and the amplifier, is charged to 
the peak value of the input signal. The amplifier output produces a con- 
tinuous voltage equal to the charge on the capacitor. Discharge time 


constants of five minutes or longer are obtainable. 


Multi-contact rotary stepping switches were used in the control 
system. Four milliseconds switching time is required with a self- 
excited stepping capability of 58 steps per second. The electrical 
noise resulting from the use of the stepping switch was not a problem 
since signal voltages were kept large whenever possible to completely 
overshacow the microvolt noise ieaqu! Also, the reader delays 


inputs until after allwipers have attained their new positions. 


A Philbrick dual power supply, Model R-100B was used to drive all 
the operational amplifiers. Typical drift for this power supply is less 
than 100 mvdc per day. 


tf 


The storage of voltages in the control system is done entirely by 
the use of capacitors. Cia and = shown in Fig. 2 are O.5uf while 
the rest of the storage capacitors are O.25uf . The charging time of 
the capacitors is ten milliseconds. After this period, the capacitors 
have attained 99.9% of their final charge. The charging time of 10Omsec 
would permit switching up to sixty steps per second. The speed in this 
control system was limited to ten steps per second by the motor drive 
on the reader head. Time base for the entire control system was gener- 


ated by an adjustable pulsing circuit shown in Fig. 4. This pulse was 


used to synchronously step the reader head and switches. 





IV. OVERALL SYSTEM ENGINEERING AND PERFORMANCE 


The four channel control system performance data is found in 
Appendix C. Tests were run to determine the accuracy and repeatability 
of multiplication by the computing system. Since accuracy of computation 
is directly effected by the peak follower output; determination of the 


R-C time constants for the storage network was also made. 


Accuracy of multiplication was found to be in error by 1.6% or less 
for computing speeds of eight steps per second. This accuracy was de- 
termined by multiplying two known reference voltages together, measuring 
the product and comparing it with the otherwise calculated value. (see 
Appendix C) 


Because of the inherent stability of the computer components and 
care taken in wiring to guard against ground loop problems, the repeata- 
bility of computation over a fifteen minute period resulted in output 


variations of less than 1.0%. 


R-C time constants for the peak follower inputs were found to be 
five minutes or greater resulting in peak follower output variations of 


less than 0.5% for a switching time of 0.12 seconds. 


Accuracy of producing cle, h , and ceon is dependent upon the 
number of discrete levels of the reader, the tolerances of the resistors 
used, and the accuracy of the reference voltage. First, 159 numbers were 
aneagleronle giving the reader a maximum inaccuracy of 0.3%. The reference 
voltage was a dry cell, and since the current demand was at most a few 
microamperes, it was considered to have constant output. Finally, the 


measured error due to resistor tolerance was +3% maximum. 


As stated previously, maximum switching speed was limited to 10 
cycles per second by the motor drive on the reader head. Since three 
steps are necessary for each determination of qs» corrected outputs 
of desired heat flow were made every 1.2 seconds for each of the four 
channels. This speed could have been increased by a factor of five 
simply by changing the motor speed without adversely effecting the ac- 


curacy of the control systen. 





Total cost of the control system described in this work was 


approximately $1,500. 


To prevent the hazards of infra red radiation exposure to operating 
personnel, the control system and instrumentation was remotely located 
in a room overlooking the test sight. <A description of the facility 
layout is given in Appendix B. 





V. CONCLUSIONS 


From the results of this investigation it is concluded that a 
time-shared digital-analog computer control system for aerodynamic 


heating is 


(1) capable of heat simulation control up to Mach 6; 
(2 sitiple. 

(3) dependable. 

(4) safe to operate. 


(5) relatively inexpensive. 
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APPENDIX A 


Analog Summation 


The operational amplifier used in the summation circuit has an input 


impedance which is assumed to be infinitely large. 





Since the input impedance is infinite, the current through terminal 


A is zero and currents through R and R 


Z ri 


equation may be written: 


Voe=x=-MV 
a 


Eliminating V, from equations (1) and (2) gives: 


Ms ue 
Re (Vv, - =) 7 -R,(¥, +) 


Or rewriting 


V RM 
oe i 
Vv.  R, +R.(1 +M) 

z a i 

If M is very large, the approximation is'made: 
R ia 
Me -- = Ye 
ae 


For several inputs this can be extended. 


mae ee 


are equal. 


The following 


(1) 


(2) 


(3) 


(3a) 


(4) 





7 
a 
>| ee 
© 





ae iL f 
ale 
Ro 
V0 
R 
5 
We 3 Wh 


R R R ae : 
vo (38) - (ri) Ye -(e) % o 


Multiplication 


ae is generated as a voltage signal by multiplying it by a reference 


voltage. 


Output of -t (Step 1) 





tis added to -t ina summer (See equation 5) to give (to7 t) 


Che oe) (step 2) 


The analogous voltage, Coes t); is stored in a capacitor and 


appears as a continuous output of a peak follower. 


é ‘ 


a t.) (Step 3) 


a 
= : 
1 
wees 
I) 
ae 


The output of the peak follower is fed back to the tape reader as a 


reference voltage where it multiplies h to give q. 


ee ee 





| Tape Reader 
h a CARE a q= (ten (Step 4) 
aw 8 
Ana log 
\ 


55, 7 t.) 
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APPENDIX B: STANFORD UNIVERSITY RADIANT HEATING STRUCTURAL TEST FACILITY 


I. INTRODUCTION 


The need for aerodynamic heat simulation in structural testing of 
supersonic aircraft is well recognized. Thus, the authors undertook 
the design and construction of a four channel radiant heating test facil- 
ity at Stanford University during the summer of 1962. The system was | 
designed to ieet the following requirements: a high intensity, rapid 
response ieat source; an adequate reflector system to direct the heat; 
an absorptive finish on the test Specimen to permit maximum utilization 
of available heating power; a practical technique for surface temperature 
determination; a system of radiation measurement; and a computer control 


system to correctly vary the heat input to the test specimen. 


II. SYSTEM LOGIC 


Figure B-l is a block diagram of the control system for the radiant 
heating test facility. Voltages analogous to h and co are generated 
and then utilized by the computer along with a voltage analogous to te 
for solution of the Newtonian convective heat transfer equation, 
gq = h(t. - t.) - Another voltage (a) representing radiation from 
external sources (sun, atomic blast, etc.) is generated and then added 
to q ina summing circuit to give the desired heat flow signal, da - 
The: comparator output, c , is the difference between dg and ae ; 
the measured heat flow. The control signal causes the power regulating 
device to vary the power supplied to the heaters. The heat flow measSure- 
ment system closes the control loop, causing the correct heat input to 
the test structure to be maintained. The measurement of heat input to 
the test specimen is made by either of two systems. The operational 
details of these devices are explained briefly in Section III and in 
detail in References [1] and [2]. 


III. SYSTEM DETAILS 


Certain requirements are essential in a practical feedback control 
system for the simulation of aerodynamic heating. The requirements set 


Hom tniS Tacility are te: 


oe eens 





(1) accurately compute desired heat flows for simulation of speeds 


up to Mach 6. 
(2) be capable of measuring the heat inputs to the test specimen. 
(3) be dependable. 
(4) be safe to operate. 
(5) be inexpensive as possible. 


(6) be simple to operate and program. 


The operations of the function generator and computer shown in 
Fig. B-1l are combined in the digital to analogue computer control system. 
This unit consists of a punched tape reader, an operational amplifier 
and resistance bank, two peak followers and associated capacitors, a 
multi-bank stepping switch, a summer, and a chopper stabilized amplifier. 
Calijeas B-2) The details of the individual components are given‘in Ref- 


erence [1] and only the modus operandi is discussed here. 


Values of h and (eae are obtained from the Stefan-Boltzmann law. 
A binary coded program of these time dependent functions is made on a 
punched tape. This tape is placed on a reader and used in the time 
Shared system as follows: an analogue voltage proportional to ae is 
generated from a resistance network and a reference voltage in conjunction 
with an operational amplifier. This signal is temporarily stored. The 
next command is t_.. which is treated in the same way. The output of 
the temperature transducers, t. » is amplified by a chopper stabilized 
amplifier and then subtracted from the stored c analogue in a Sunmer. 
The voltage difference COR = t.) is used as the reference voltage for 
the next stage of the operation. This consists of giving the command, 
read h , and by virtue of the fact that the reference voltage is 
analogous to cone t.) the output of the system is analogous to the 
product h(t - t.) - Simultaneously, q and qa, are added in a 
summer and this voltage, analogous to qq? is Stored in a capacitor. 
These steps are repeated for channels two, three, and four. The digital 
to analog computer evaluates a new value of da for each channel every 
1.2 seconds. ‘Thus, by time sharing, by digital to analog conversion, 
and by unique computer Ceemmauee Oe: a Simple and accurate means for 


obtaining a da history is realized. 


Sees 





There are two systems available for closing the control loop. Both 
systems meaSure the rate of heat flow to the test specimen and therefore 
may be used interchangeably. One system, the classic system, consists of 
special transducers. These essentially consist of thermocouples in con- 
junction with substantial heat sinks. Ina radiant field the thermo- 
couples reach an equilibrium temperature proportional to the radiant 
intensity. Each transducer is two thermocouples connected so that their 
signals oppose. One is subjected to the radiation from the heaters and 
the other to the reflected radiation. The output is therefore proportional 
to the difference between the incident and reflected radiation. This is 
qa exclusive of convective losses. The heat sinks for the transducers 
are cooled metallic blocks. The instruments are attached directly to the 


Sides of the heater arrays. They are calibrated in a calorimetric manner. 


_ The other system, one especially devised for this project, consists 
of eighteen thermocouples arranged to obtain temperature profiles in the 
test specimen. <A high speed stepping switch sequentially connects the 
thermocouple signals to an analogue computer which solves the appropriate 
heat flow equation. The computer is time shared by the 4 stations. The 
latter system is more versatile than the former since it measures net 
heat flow irrespective of the nature of the source. Thus convective losses 
in Simulation are accounted for. An additional refinement to the computer 
permits compensation for radiant losses which the structure would normally 
experience in flight. Figure B-3 is a photograph of the feedback com- 


puter system. 


Voltages in the control system are scaled such that at any time the 
prescribed qs level is twice the anticipated q, level. A summing 
circuit subtracts q from q, and the output (c) is used to control 
the power regulating devices. Due to the two-to-one scaling the control 
signal consists of a set-point voltage plus an error voltage. Thus, the 
heaters are controlled in such a manner that the desired and actual heat 
inputs are equalized. This occurs when the error voltage level converges 


to zero and the control signal converges to the set-point voltage. 


Main power is supplied through two 300 KVA, 480 volt, 3 phase, 
General Electric transformers. The power regulating devices consist of 


75 KVA, 480 volt, 3 phase General Electric saturable core reactors and 


ee 





their associated control units. Transient operation of the power equip- 
ment on an overload duty cycle makes possible peak outputs Gr two tO 
three times rated value. Heating capacity of more than one megawatt is 
therefore feasible. Figure B-4+ is a photo of the saturable core reactor 


nits . 


The basic heating elements are General Electric 7.5 infra-red lamps. 
(See Fig. B-5 ) They consist of a tungsten filament supported by tanta- 
iim disc spacers in a sealed quartz tube ey of an inch in diameter and 
12 inches in length. Each unit is rated at 1000 watts at 220 volts. 

They are capable of withstanding the 480 volt peak output of the power 


supply and under these conditions develop approximately 3 kilowatts. 


The heaters are arranged in groups of six. Figure B-6 isa photo- 
graph of one such unit. The heater arrays are connected in delta to the 
three phase power from the saturable core reactors. Each unit covers 
about one third of a square foot. Thus, the maximum radiant intensity 
possible with the present arrangement is about 60 kilowatts per square 
“foot. Two thousand watt heaters of the same external size are avail-' 
able when higher intensities are required. Peak output of these is 6000 
watts per heater. ‘The heater arrays are quite adaptable and may be 
bolted to any suitable frame to surround a test structure. Special 
arrangements may easily be devised for special purposes, such as the 


heating of curved surfaces. 


Figure B-7 shows an overall view of the test facility. The main 
items of equipment and the channels of information flow are outlined on 
the photograph. All control, computing and recording equipment is located 
in the control room. This room is on the second floor of the laboratory 
Overlooking the test site. Large windows provide for visual observation 
of the structure under test. The remote location of the control room 
was dictated by safety considerations and provides a convenient, compact 


Center for operation of the entire facility. 


IV. CONCLUDING REMARKS : 


The design, construction, and operation of this test facility has 
proved that the concept of component time-sharing is practical in certain 


types of control systems. It is concluded that such time-sharing is an 
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economical proposition in which savings in equipment costs are balanced 
against slightly reduced system performance. The facility is such that 
it can be readily installed in any small research laboratory where cost 


is of prime importance. 
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Fig. B-4 
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Fig. B-5 





Heater Array 
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Reference Voltage Multiplier Actual Prod. % Error 
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1.88 volts 3.60 6.70 volts 1610 
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Seca VOLES 2.20 Sess) areas 1-2 
1On@ sn volts 2.20 42.1 volts 0.5 
Duplication of Above Test after 15 Minutes: 
iS Onvokrs 3.60 Be . aigoyll ae 0.9 
1.88 volts 3.60 Gane volnes 120 
2.02 volts 3.10 6.28 volts OF) 
Sr lO ve las Saale 9.61 volts 0.0 
4.30 volts 3.00 12.38) volts 0.8 
more® volte 3.00 18.4 volts Lae 
Oe veils 2.95 24.5 volts 156 
10.4 volts 2.95 SOS arelhes 0.6 
12.4 volts 2.60 Sl OGmevoltuc lie 
ae) ideas: 2.60 35.5 volts Onc 
(22 VOCS 2.20 Soagm Volts 0.3 
19.0 volts 2.20 42.5 volts gore 
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